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Abstract 

We have developed a Monte-Carlo simulation code for an aerogel Cerenkov 
Counter which is operated under a strong magnetic field such as 1.5T. This code con- 
sists of two parts: photon transportation inside aerogel tiles, and one-dimensional 

amplification in a finc-mcsh photomultiplicr tube. It simulates the output photo- 
electron yields as accurately as 5% with only a single free parameter. This code is 
applied to simulations for a B-Factory particle-identification system. 

PACS number : 29.40.Ka 

keywords : aerogel, Fine-Mesh photo multiplier tube, Monte-Carlo 
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1 Introduction 



Particle identification (PID), particularly the identification of charged pions and kaons, 
plays an important role in CP-violation studies in B-Factory experiments In the 
BELLE experiment at the KEK B-factory(KEKB), a threshold silica aerogel Cerenkov 
counter (ACC) will be used to extend the momentum region of the PID capability beyond 
the reach of dE/dx and time-of- flight (TOP) measurements]^]. 

Silica aerogels (aerogels) are a colloidal form of glass, in which globules of silica are 
connected in three-dimensional networks with siloxan bonds P, Many high-energy and 
nuclear-physics experiments have used aerogels instead of pressurized gas for Cerenkov 
counters. 

In the BELLE experiment, fine-mesh photomultiplier tubes (FM-PMTs) will be used 
for the readout of Cerenkov photons emitted from aerogel]^, |^. This was only choice at 
the design period for this kind of counter operated under a strong magnetic field, such as 
L5T. 

We have developed a Monte-Carlo simulation for the aerogel counter, which includes 
photon transportation in the aerogel tiles and response of a phototube. One of the key 
issue in this simulation is the parameterization of absorption length, which is obtained 
by a careful comparison of beam test data and the prediction from the simulation. We 
describe a Monte-Carlo simulation of this detector system in this paper. In the next 
section we introduce our experimental setup of the threshold Cerenkov counter. The 
purpose of this Monte-Carlo project is also described. The framework of this Monte- 
Carlo code is described in the third section. The simulation of fine-mesh phototubes is 
described in section 4, and photon transportation inside aerogel tiles is shown in section 
5. In section 6, a simple model for aerogel absorption is introduced. The results on this 
model is described in section 7. In section 8, we discuss on the wavelength dependence of 
the light absorption. The conclusion is given in section 9. 
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2 Threshold Cerenkov Counter with a Fine-mesh Pho- 
tomultiplier Tube Readout 

We will use threshold Cerenkov counters to identify charged kaons in the momentum 
region 1-4 GeV. The configuration of the aerogel counter system of the BELLE detector 
is shown in Figure ^ It should be noted that the KEK B-factory is an asymmetric 
e^e~ collider with energies of 3.5 x 8 GeV and the produced particles from the B meson 
decay are Lorentz boosted to the e~-beam direction(z). Since the momentum range of 
the produced pions and kaons changes as the emitted polar angle changes in a laboratory 
system, we employ the different refractive indices (n) of aerogel tiles corresponding to the 
different polar angle. In the endcap region, however, n=1.03 was selected for i?-flavor 
tagging purpose 0. 

The typical arrangement of the counter is shown in Figure H. The size of an aerogel 
tile is approximately 12 x 12 x 2.4 cm^. The single counter contains 4, 5, or 6 layers of 
them with a total thickness of 10-14.5 cm. The Cerenkov photons are readout by one or 
two fine-mesh photomultiplier tubes (FM-PMT). In order to reduce effects of a magnetic 
field of 1.5T in the z direction, we selected the FM-PMT. We will use three types of FM- 
PMTs of 2, 2.5, and 3 inches in diameters. The effective diameters of these FM-PMTs 
are 39, 51, and 64 mm. The inner side of the counter is covered with a white reflector 
(GoretexH). 

The dominant scattering mechanism of the Cerenkov photon inside the aerogel tile 
is considered to be Rayleigh scattering. There is almost no understanding about the 
absorption mechanism of the light inside aerogel. The main purpose of developing this 
Monte-Carlo code is to understand this mechanism, or in other words, to construct a 
reliable model. 

The ampliflcation mechanism of the FM-PMT in a strong magnetic field is also poorly 
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understood. Since the B-factory is a high-statistics experiment, an accurate Monte-Carlo 
simulation is necessary for these PMT's responses. 

In summary, we would like to develop a Monte-Carlo simulation with an accuracy of a 
few percent level. It should simulate such effects as position dependence and edge effect. 

3 Framework 

The typical scattering length of aerogel is a few cm at A =400nm Therefore, Cerenkov 
photons undergo many scatterings (including absorptions and reflections) before reaching 
photocathodes. It is very difficult to understand this situation analytically. The Monte- 
Carlo simulation may be the only calculation choice. 

In coding the Monte-Carlo simulation, we used C++ language. A part of GEANT3.210 
was used for the geometry and display parts. 

The flow chart of this Monte-Carlo simulation (hereafter, referred as "ACC++") is 
shown in Figure ^ Each box corresponds to a major process in the program, which 
is written as "class" in C++. The vertical arrows indicate the direction of the data 
flow. The display and geometry parts are common utilities. At first the hit points by 
the charged particles inside the aerogels are recorded. Cerenkov photons are emitted at 
the above-mentioned hit points according to the velocity of the charged particles and 
the refractive index of medium. Then, the photon transport part takes care of photon 
transportation to the phototube window surface. Finally, the FM-PMT part amplifies 
the photoelectrons and the ADC counts are obtained. For data management, we used the 
CLHEP library 0. 

4 Fine-mesh Photomultiplier-Tube Simulation 
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4.1 Structure of a Fine-mesh PMT 

A sectional view of the FM-PMT is shown in Figure ^ Each PMT has a borosihcate glass 
window, a bialkari photocathode, 19 fine-mesh dynodes and an anode. The diameters of 
the effective area(0), equal to the diameters of the dynodes, are 39, 51, and 64mm for 
2", 2.5", and 3" PMTs, respectively. The cathode-to-anode distances(L) are 20, 20, and 
23mm for 2", 2.5", and 3" PMTs, respectively. The FM-PMT used for the BELLE ACC 
has the following improved performance. The average quantum efficiency(QE) of the 
photocathode is 25% at 400nm wavelength, improved for recent products |ll|. A finer 
mesh than that of conventional products is used to improve the gains in the magnetic 
fields 0. The optical opening of the mesh is about 50% [Q. 

4.2 Simple Model of Amplification 

The first step of the simulation is conversion of incident photons to photoelectrons(p.e.). 
This conversion is simulated while taking into account the wavelength dependence of 
the quantum efficiency (Qi?). Numerical data of the wavelength dependence, obtained 
from the manufacturer, are used. The absolute value of QE for each PMT is deter- 
mined by the cathode blue sensitivity, denoted as S'A;6[|lT|. The absolute QE at 400nm 
wavelength(Qii^(400nm)) can be calculated as 

QE(400nm) = 0.026 x Skb (1) 

using the data obtained from the manufacturer. The QE does not depend on the diameter 
of the PMT. The average Skb value is 9.7. 

The second step is to simulate the amplification in fine-mesh dynodes. A simple model 
is schematically shown in Figure ^. For each dynode stage, some incident electrons pass 
through holes of a fine-mesh dynode, and others hit the dynode wires, followed by the 
emission of secondary electrons. The probability of passing through holes is proportional 
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to the optical opening(the ratio of the hole area to the total area) of the mesh, /(~ 0.5). 
Some of the secondary electrons pass through the holes, and others are re-absorbed by 
the same wire. The re-absorption probability is dependent on the hit position, and it is 
higher when electrons hit the mesh at the top of wires, compared to the case when they 
hit it at the edge. That probability increases at higher magnetic fields. This dependence 
is considered in the simulation by introducing a blind region, where all of the secondary 
electrons are assumed to be re-absorbed. The ratio of the blind region area to the total 
dynode area is represented by a parameter b in the simulation. The effective gain per one 
stage of a djnode{6^^^) is then parameterized as 

S'^ff = f + 6'^'\l~f-b), (2) 

where 5^'* is the number of secondary electrons passing through holes, whose distribution 
is assumed to be a Poisson in the simulation. The quantity 5''** is determined so that the 
total gain for the 19 stages(G') is reproduced. Since the gains of the first three and the 
last two dynodes of our FM-PMTs are improved by a factor of A;||l2|, we obtain 



Here, S'^^'^ is typically around 2.5 . 



In a simulation of the pulse height for each incident photoelectron, the above process 
is repeated until the sixth dynode. It is shown later that a calculation up to the sixth 
dynode is sufficient to reproduce the pulse-height spectrum. 

4.3 Single-photoelectron Spectrum 

Using the method described in section 4.2, a pulse-height spectrum for single-photoelectron 
events is simulated, as shown in Figure (a). That is compared to the experimental data 
shown in Figure ^(b). The experimental data were taken by using a pulsed laser|13]. The 



number of events in the pedestal region is more than 90% of the total events. Therefore, 
most of the signal events are single-photoelectron events. 

The simulation reproduces the experimental data well. The single-photoelectron spec- 
trum of FM-PMTs has no characteristic peak, in contrast with that of line-focus PMTs. 
Only a fractal structure is observed in the spectrum, the interpretation of which is de- 
scribed in reference ^j. 

4.4 Multi-photoelectron spectrum 

As noted in section 4.3, a FM-PMT has no peak corresponding to a single photoelectron 
in the spectrum. This deteriorates the pulse-height resolution of a multi-photoelectron 
spectrum, compared to the resolution determined by the original Poisson distribution of 
incident photoelectrons. 

Figure |^ shows a comparison of the incident photoelectron distribution(a) and the 
simulated spectrum at the anode for corresponding events(b). To quantitate the deterio- 
ration, we define the excess noise factor(ii^A^F) as 



where fiinput is the average number of input photoelectrons and /^e// is the effective number 
of output photoelectrons. In the present work, the output spectrum is fit to an asymmetric 
Gaussian with two width parameters (ai and (J2), each representing the width of the lower 
and higher sides, respectively, as shown in Figure 0-(b). Then /Zg// is calculated as 



where n is the mean of the spectrum (mean ADC count). In the example shown in Figure 
0, fJ'input is 20.0 and /Xe// is 9.64. Therefore, ENF is 2.07. 



ENF 



(4) 




(5) 
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Figure | shows the measured ENF for 2", 2.5", and 3" PMTs. The tested PMTs 
were illuminated by the same 412nm pulsed laser light. The number of incident pho- 
tons was calibrated using a reference PMT, that is a 2" line-focus type PMT (Hama- 
matsu R329-05S), which was well calibrated by the manufacturer. The average num- 
ber of photoelectrons(yuf^*" J was deduced from the average ADC channel(g) and gain. 
The gain was obtained from the mean ADC channel for a single-photoelectron event(go)- 
Taking into account of the difference in the quantum efficiency between the reference 
PMT(QE(ref)) and the tested FM-PMT(QE(FM)), the number /if^*^^ was obtained as 

^ QE{¥M) q 
^^input QE(ref) go' 

The quantum efficiencies {QE{FM) and Qi?(ref)), measured by the manufacturer have 



about a 10% ambiguity, giving about a 14% error in the measured i?A^F|]T2|. The mea- 
sured ENF is 1.99 on the average. The agreement between the simulation and the data 
is satisfactory. 



4.5 Effect of a Magnetic Field on the Resolution 

The electron trajectories are influenced by magnetic fields. In the simulation, the effects 
of magnetic fields are treated as follows. 

On the microscopic scale, electrons spiral along the magnetic field. At 1.5T, the Lamor 
radius of electrons for a typical energy of 6eV is 5.5/im, which is comparable to the mesh 
width. Therefore, the probability for the re-absorption of secondary electrons by mesh 
wires increases. This effect can be regarded as an increase in the blind region(6) on the 
mesh, which is shown in Figure ^, not only as a decrease in 5^***. The increase in b affects 
the pulse-height resolution, while a decrease in 5^** affects the total gain. 

On the macroscopic scale, electrons move parallel to the direction of the magnetic 
field. Therefore, some electrons can not reach the anode when the PMT axis is inclined 
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with respect to the magnetic field. The fraction of such electrons is determined by the 
diameter of the photocathode(0) and the cathode-to-anode distance(L). In the simulation, 
the position of secondary electrons on the anode plane is calculated from information 
about the positions of electrons at the cathode and the field direction. Then, in case that 
an electron is in the effective area of the anode, the amplification process is simulated. 

Table [l| shows the ratio of the effective number of photoelectrons with magnetic fields 
of l.STesla (/ie//(1.5T)) to that without magnetic fields (/ie//(OT)), obtained from our 
simulation and a measurement for each PMT size. The results are shown for 6 = 0° and 
35°, where 6 is the angle between the PMT axis and the fields. The parameter b is adjusted 
to reproduce the experimental data for each angle. In the case of ^ = 0°, the experimental 
data can be reproduced well with b = 0.15. In the case of 6 = 35°, the results of the 
simulation are shown for two cases, 6 = and b = 0.15. A better agreement between 
the data and the simulation is obtained when only the macroscopic effect is considered 
(6 = 0). Therefore, it is indicated that the macroscopic effect dominates for 6 = 35° case. 
The dependence of the parameter b on the angle 6 can be explained by a change of the 
mesh width projected along the magnetic field direction. 



5 Photon Transport inside Aerogel Tiles 
5.1 Algorithm 

A charged particle is transported inside aerogels by a 1-mm step. Cerenkov photons are 
radiated along the charged-particle's trajectory. The number of photons (Np) is given by 
the Frank- Tamm equation [Q, 
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where a is the fine-structure constant, L the thickness of the radiator, z the charge of 
an incident particle, 6c the Cerenkov angle, and E the energy of a radiated photon. The 
Cerenkov angle (6^) for a particle with velocity f3c in a medium with index of refraction 
(n) is 

Oc = cos-^(-l). (8) 

pn 

A radiated photon is transported according to the algorithm shown in Figure |^. In 
the photon-transport part, we first obtain information about photon, i.e., the wavelength, 
radiated position, and radiated direction from the Cerenkov-radiation part. 

The Cerenkov photon is transported to the next position in the following way. The 
mean transport length (A) is defined as 

111 

+ ^^^ (9) 



A(A) AscatW AabsW 

where A^cat and Aabs are the scattering length, and absorption length respectively, in 
the aerogel. They are calculated based on the spectrophotometer measurement, which is 
described in section 5.2.1. The transport length is calculated using the exponential prob- 
ability function (Poisson distribution). It is then judged that the position after transport 
is (i) outside the counter box or (ii) inside with the information from the Geometry part, 
in which we used GEANT 3.21 (only geometry routines) 0. (i) If the position is out of the 
box, the photon is traced to the inner wall of the box. We then judge whether the position 
is on the FM-PMT surface or diffuse refiector surface (Goretex 0). If the position is on 
the FM-PMT surface, the photon information (the wavelength and the position of hit) is 
given to FM-PMT simulation part and the photon transport is finished. If the photon 
hits a refiector surface, we judge whether the photon is absorbed or refiected by Goretex 
based on a spectrophotometer measurement of the Goretex refiectivity, which is described 
in section 5.2.2. If the photon is absorbed by Goretex, the photon transport is finished. 
If the photon is refiected by Goretex (diffuse reflection), the photon is transported again. 
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(ii) If the position after transportation is in the counter box, it is judged that the photon 
is absorbed or scattered by aerogel, according to the following ratios: 

-^'■scat ~r ^^abs 

Pscat = T J^T , (11) 
'■scat I J- i-abs 

where the Pats and Pscat are the probability of absorption and the scattering of the aerogel, 
respectively. If the photon is absorbed by aerogel, the photon transport is finished. If 
the photon is scattered by aerogel, we randomize the direction according to the Rayleigh- 
scattering formula and repeat the above-mentioned transportation again. 

5.2 Input Parameters 

5.2.1 Spectrophotometer Measurement 

We use the spectrophotometer [jl5| in order to measure the transmittance of aerogel and 



reflectivity of Goretex as a function of the wavelength(250 ~ 800 nm). The transmittance 
(T) is defined as I/Iq, where / and Jq are the photon fiuxes at the photodetector with 
and without aerogels, respectively. The thickness of a typical aerogel tile is 24mm. The 
spectrophotometer has a finite solid angle at the entrance of the photo detector. A part 
of scattered photons is detected by the photodetector. Therefore, the transmittance (T) 
is 

T^l-A, (12) 
Ty^l-A-S, (13) 

where A is the absorption and S is the scattering. In order to estimate the contribution 
from S and A to T, we took the geometry of the photospectrometer and other necessary 
information into account for our Monte-Carlo (ACC++) and carried out a simulation. 
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The details are discussed in section 5.2.3. A simulation display of the photospectrometer 
is shown in Figure |10[ The lines are the photon trajectories. The circle is the aperture of 
the photodetector (8mm0). 



5.2.2 Goretex Reflectivity 

It is not possible to measure the absolute reflectivity of a material directly with the spec- 
trophotometer. We therefore measure the relative reflectivities of samples to a reference 
sample whose absolute reflectivity is known. We used an NIST traceable standard reflector 



(Spectralon) which was well calibrated by a company|16|. Then, the absolute reflectiv- 
ity of Goretex was calculated. Figure |ll| shows the measured absolute reflectivity of the 
Goretex as a function of the wavelength. The reflectivity of Goretex is better than 93% 
over the range of the measurement. It shows a better reflectivity in the short- wavelength 
region, i.e., a better acceptance for Cerenkov hght (UV light). 

5.2.3 Rayleigh Scattering 

We input two parameters [Aabs and Agcat) in the simulator. Initially, we suppose that the 
transmittance is only due to the Rayleigh-scattering effect. We consider the absorption 
and scattering parameter, as expressed below by Equations |T^ and |T3|, and input them to 
the simulator: 

Aabs = OO, (14) 

—d 

Ascat = (15) 

where d is the thickness of the samples. Although the simulated transmittance shows 
almost the same shape as the experimental transmittance, the absolute value of the sim- 
ulated data points are higher than the experimental data (i.e., more transparent). This 
indicates that the scattered events within the finite solid angle of the photodetector and 
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the absorption events are taken into account for the transmittance. We must take these 
factors into account. Then, the aperture factor (/„) is introduced in Equation |l^ and the 
definition of ^scat is modified to 

Kcat = faX (16) 

To evaluate the aperture factor {fa), a simulation was carried out. The obtained value 



for fa is 0.55. Figure shows the transmittance data and simulation result, which agree 
with each other. In other words, the photospectrometer measurement of our aerogels 
suggests that there is no, or very small, absorption. 



6 Model for Absorption 

Empirically, it was known that the transmittance obtained by the spectrophotometer was 
well fitted with the two terms: 

Ascat = aX'^, and Aa^^ = bX'^, (17) 

where A is the wavelength and a and b are free parameters. The term was empir- 
ically introduced in order to improve the fitting of the transmittance spectra. In the 
reference ||17||, they assumed that the A'^ term is the Rayleigh scattering, and that the A^ 



term is the absorption. The results (n=1.010) are shown in Figure |T3|for our aerogel. 

The typical scattering length was calculated to be ~ 2.5 cm at 400 nm wavelength, 
and the absorption length was fitted to be 10-times the scattering length. When we input 
parameters Agcat and Aabs into our Monte-Carlo simulation (ACC++), we obtained only 
~ 2 photo-electrons (p.e). On the other hand, an experiment using 3.5 GeV vr" beam was 
carried out for our counter box (Figure |^); we obtained ~20 p.e., which is inconsistent 
with the above assumption. Therefore, we conclude that bX^ is not due to absorption, 
and that the absorption length is greater than 100-times the scattering length. 
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In order to estimate the magnitude of the absorption length, we make a simple model 
that the absorption length is proportional to the scattering length, i.e., we assume the 
same wavelength dependence as that of the Rayleigh scattering. 



Kbs = a X A 



scati 



where a is a free parameter. The a value was determined by comparing the experimental 
data and simulation results as described in the next section. We then input these two pa- 
rameters {Aabs and Agcat) into ACC++. Figure shows a simulation display of ACC++ 
for the photon transportation inside the aerogel tile with the assumption a=100. The 



other type of absorption model can be found, for example, in reference|18 . 



7 Comparison between the Beam- Test Results and 
Simulation 

7.1 Determination of a by Minimization 

We carried out a minimization by comparing the fieff value of the experimental data 
and that of the Monte-Carlo data while changing the free parameter {a). The experi- 
mental data were obtained in the beam test with the 3.5-GeV tt^ at the it2 beam line 



of the KEK PS. The incident position of the vr 's is the center of the box. Figure 15 



shows this result. From the minimum of this curve we obtained a = 408.1±5g f which 
corresponds to an absorption length (Aabs) of about 5.4 m at 400 nm. The a values differ 
(by a factor of ~2) for different samples of aerogels having the same refractive index and a 
different refractive index. There are no clear dependence on transmittance nor refractive 
indices. We, therefore, must obtain this value for each counter box, for example using 
e~^e~ —>■ fJ^~^fi~ events in the real B-factory experiment. 
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7.2 Position Dependence and Accuracy 

With the best-fitted a value we compared the beam data and the Monte-Carlo prediction 
for the position dependence of light yields. The experimental data of /ig// were obtained 
for the 3.5-GeV vr^. The data were taken at various positions of the incident beam. 
We selected 23 incident positions. The coordinate system of the beam incident position 



and direction are defined in Figure 16. The beam direction is (0, 0, 1) for all data. A 
comparison with the Monte-Carlo data are shown in Figures 0: (a) at y=0 plane, (b) 
y=2cm, (c) y=3.5cm, and (d) y=4cm. The hatched areas are the Monte-Carlo predictions. 
These errors are statistical. The disagreements between the experiment and Monte-Carlo 
values are around the 5% level. We therefore consider the accuracy of our Monte-carlo 
simulation to be 5%. 



7.3 Miscellaneous Distributions 

With the best-fitted value of the parameter a, we derived the distributions of the number 
of scatterings inside the aerogel which our Monte-Carlo simulation predicts. They are 



shown in Figure |1^. The mean number of reflections on the Goretex is ~ 11-times, that 
of the scattering in aerogel is ~ 29-times, and mean total transport length is ~ 83 cm. 
43% of initial radiated photons in aerogel reach the PMT surface, 32% are absorbed by 
Goretex, and 25% are absorbed by aerogel. This Monte-Carlo code is very useful in 
designing the shape of the counter box for future experiments. 



8 Wavelength Dependence of Absorption 

As noted in section 7.2, our Monte-Carlo simulator is able to reconstruct beam data with 
5% accuracy after fitting the a value. For a further study to evaluate the wavelength de- 
pendence of the absorption length(Aab5), we carried out experiments with several types of 
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optical filters (sharp-cut filters ||T9[). These filters cut the light with a shorter wavelength 
than threshold values. Here, we define the filter's name, whose threshold wavelength is 
300nm, 300nm-filter. We used eight types of filters, whose threshold wavelengths are 
300nm, 360nm, 400nm, 460nm, 500nm, 560nm, 600nm, and 700nm. The transmittance 
of each filter was measured with spectrophotometer as a function of the wavelength and 
input to ACC++. 

The beam-test setup and evaluation procedure of the a value were the same as those 
mentioned in section 7.1. The refractive index of aerogel was n = 1.013. We used a 3" 
FM-PMT with a 2" mask, because the diameter of the filters were 2". The filter was 
set between the PMT surface and the aerogel tiles. After fitting beam data with the 
simulation results, we obtained the best a value for each filter setup, a = 302.7±gg for 
the 300nm filter, a = 596.8±m for the 360nm filter, a = 3104.5±^?J:^ for the 400nm 



filter and a diverges for the longer-wavelength filters. These are plotted in Figure 19 



function of the threshold wavelength of the filter. From these a values, we can conclude 
that Aabs is not proportional to A^. The absorption is significant only in short-wavelength 
range (around 300nm), but not in the longer wavelength region (>400nm). 

Another absorption model, which assumes Aabs = exp(/5A), is found in the reference 
T8| . Here /3 is a free parameter. An attempt to evaluate (3 was carried out using the 



minimization with the center incidence data. Then we looked at the wavelength depen- 
dence (experimental data with the sharp-cut filters) of the Np^. In the result, we find that 
f3 is still dependent on the threshold wavelength of filter. The (3 value increases at longer 
wavelength. This again suggests that the absorption only happens at short wavelength 
such as 300nm. 
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9 Conclusion 



We have developed a Monte-Carlo Simulation for an aerogel Cerenkov counter under a 
strong magnetic field. This code consists of two parts: photoelectron amplification in a 
fine-mesh phototube, and photon transportation inside aerogel tiles. 

In the former, we have developed a simple one-dimensional fractal-amplification model. 
It successfully reproduces the deterioration of the resolution and the single-photoelectron 
response in a fine-mesh phototube. 

In the latter, we have proven that recently produced aerogels have a significantly longer 
absorption length, such as longer than a few meters. Here, a simple model agrees with 
the experimental observations. This model only uses a single free parameter. 

With this Monte-Carlo simulation code, we can simulate effective photoelectron yields 
within 5% ambiguities for various shape of counters, whose geometry can be defined in 
the GEANT framework. This program will be used for a particle-identification system 
for the KEK B- Factory, BELLE experiment. 
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Table R.Suda et al., NIM-A 



Table 1: Ratio of the effective number of photoelectrons with magnetic fields of 1.5 Tesla 
(/ie//(1.5T)) to that without magnetic fields (/ie//(OT)) for 2", 2.5", and 3" PMTs for our 
simulation(MC) and experimental data(Data). The detail on /ig// can be found in the 
text. 
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Figure @, R.Suda et al., NIM-A 

Fine-mesh PMT 




Figure 2: Typical form of the aerogel Cerenkov counter. 
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Figure 0, R.Suda et al., NIM-A 




Figure 3: Flow chart of the Monte-Carlo simulation code (ACC++). 
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Figure 3, R.Suda et al., NIM-A 
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Figure 4: Sectional view of FM-PMT. 
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Figure 0, R.Suda et al., NIM-A 
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Figure @, R.Suda et al., NIM-A 




Figure 6: Single photoelectron spectra: (a) simulated spectrum and (b) experimental 
data. 
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Figure 0, R.Suda et al., NIM-A 




25 30 35 
ADC(arbitrary unit) 



Figure 7: Simulated multi-photoelectron spectra: (a) the spectrum of incident photoelec- 
trons and (b) output photoelectrons. 
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Figure @, R.Suda et al., NIM-A 
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Figure 8: Measured Excess Noise Factor (ENF) for 2", 2.5", and 3" PMT. 
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Figure 0, R.Suda et al., NIM-A 
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Figure 9: Algorithm of the photon transport part of the Monte-Carlo (ACC++) 
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Figure [tO], R.Suda et al., NIM-A 




Figure 10: Event display of the spectrophotometer simulation. The lines are the photon 
trajectories. The circle is the aperture of the photodetector. 
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Figure [tl], R.Suda et al., NIM-A 
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Figure 11: Reflectivity of tlie Goretex as a function of the wavelength obtained by the 
spectrophotometer. 



32 



Figure [t2], R.Suda et al., NIM-A 
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Figure 12: Data for the transmittance of n = 1.010 aerogel (2.4cm thick). The sohd hne 
is the simulation result with only Rayleigh scattering. 
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Figure [t3], R.Suda et al., NIM-A 




Figure 13: Transmittance(circles) measured by the spectrophotometer for n=1.010, 2cm- 



thick aerogel; the sohd curve is the best fit with the parameterization (Equation |T^, the 
dotted hne contribution of absorption, and the dashed hne that of scattering. 



34 



Figure [t4|, R.Suda et al., NIM-A 




Figure 14: Event display of the Monte-Carlo simulation. The lines are the photon trajec- 
tories. 
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Figure [tS], R.Suda et al., NIM-A 




Figure 15: versus a. The details of parameter a can be found in the text. 
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Figure [tg, R.Suda et al., NIM-A 
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Figure 16: Definition of tlie incident beam position. Tlie beam direction is parallel to tfie 
z-axis. 
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Figure [t7|, R.Suda et al., NIM-A 
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Figure 17: Position dependence of /ig// as a function of the incident beam position: 
y=Ocm, (b) y=2cm, (c) y=3.5cm, and (d) y=4cm. The points with error bars are 
experimental data and the hatched areas are the Monte-Carlo predictions. 
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Figure [tS], R.Suda et al., NIM-A 
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Figure 18: Various distributions predicted by the Monte-Carlo simulation: a) distribution 
of the number of reflections on Goretex; b) number of scatterings in aerogel; c) total 
transport length inside the counter box; and d) wavelength distribution of initial photon 
(single hatch), photon after transport to FM-PMT(cross hatch), and photon after cut of 
FM-PMT's quantum efficiency (dense area). 
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Figure [tg, R.Suda et al., NIM-A 
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Figure 19: Evaluation of the factor a using sharp-cut fihers. The horizontal axis is the 
threshold wavelength of the filter. The vertical axis is the best a values, which were 
obtained by comparing the Monte-Carlo simulation data with the beam data. 
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